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Qh! Haba and Kleinert to non-autonomous mechanical systems, in the 

case of the time-dependent harmonic oscillator. In comparison with 
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1 Introduction 

The fundamental goal of stochastic quantization PQ-JS] is to reproduce the 
whole content of quantum mechanics by means of classical equations sub- 
jected to random perturbations. Presently, there remains the controversy 
about the physical origin of the noise, which is supposed to arise, for instance, 
from a fluctuating metrics or from fluctuations of the vacuum electromag- 
netic field. In spite of this, stochastic quantization is an attractive alternative 
for all physicists felling not so comfortable with the epistemological content 
of the traditional interpretation of quantum mechanics. 

Recently, Haba and Kleinert [U] have proposed a new approach for sto- 
chastic quantization, hereafter referred to as the HK method. HK is based 
directly on the use of Newton equations with the presence of noise. After in- 
troducing an auxiliary field, defined in terms of the solutions for the stochas- 
tic dynamical equations, a deterministic partial differential equation is found 
averaging over all stochastic processes. This deterministic equation, after a 
suitable transformation, yields the Schodinger equation for the mechanical 
system under consideration. Solving the Schrodinger equation reproduces 
the spectrum of the system. In comparison with other stochastic quanti- 
zation methods like Nelson's stochastic mechanics [7j or the Parisi and Wu 
method [HI E] , the HK approach seems to be remarkably concise. Moreover, 
as shown in jH] , the HK method can be applied to arbitrary time- independent 
one-dimensional mechanical systems. 

In view of the elegance and conciseness of the HK method, it is valuable 
to extend it to more general dynamical systems. Several avenues are open 
in this regard. For instance, one can address the generalization of HK to 
higher dimensionality, to composite systems, to systems exposed to external 
electromagnetic fields and so on. Here we consider the application of HK to 
a non-autonomous system, namely, the time-dependent harmonic oscillator 
(TDHO), characterized by a frequency function variable along time ^U]. We 
refrain from listing the large list of applications of the TDHO, including such 
several fields as quantum optics, cosmology, non-linear elasticity and hydro- 
dynamics [TT] . Here we limit ourselves to apply HK to this class of systems, 
showing that non-autonomous systems are also amenable to stochastic quan- 
tization via the HK method. We will see that the time-dependence of the 
frequency can be eliminated through a rescaling transformation. 

The paper is organized as follows. In Section II, we apply the HK method 
to the TDHO equation. In this Section, we both present the HK method with 



some more details than in the original reference jB] and apply it, obtaining a 
stochastic differential equation simulating the TDHO. In Section III we ob- 
tain the formal solution for this stochastic differential equation and in Section 
IV this formal solution is used to describe the time-evolution of the mother 
field, one of the fundamental elements in the HK method. We then derive 
the Schrodinger equation associated to the TDHO. Section V is devoted to 
the conclusions. 

2 Stochastic Differential Equation 

Our purpose is to apply the HK method to obtain the quantization of the 
TDHO equation, 

x + co 2 (t)x = 0, (1) 

where u(t) is the time-dependent frequency. The HK method consists of 
three steps, (a) postulate of a classical, stochastic equation, which reduces to 
the Newton equation of the system in the limit of zero noise; (b) introduction 
of an auxiliary field (the mother field) defined in terms of the solutions of the 
original, stochastic equation; (c) averaging over stochastic processes, yielding 
a deterministic equation equivalent to the Schrodinger equation. For the 
TDHO, we choose the stochastic equation 
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f(t)n, (2) 



where r = (x,y,0) and ft = (0,0,0), for constant Q, and where n = 
(1, 1, 0)/\/2 is an unitary vector in the diagonal direction in the plane. The 
function p is any particular solution for the Pinney's ^21 equation, 

Q 2 
p + u 2 (t)p = —, (3) 

where u>(t) is the time-dependent frequency characterizing the TDHO. It is 
convenient to take Q ^ so that p never vanish. Finally, f(t) is a stochastic 
variable with zero mean over statistical processes, 

(/(*)) = 0, (4) 

and with correlation function 

(f(t)f(t')) = h5(t-t'), (5) 



so that there is a white noise stochastic process. 

The choice for the classical dynamical equation is justified because, treat- 
ing formally / as an ordinary function, the second order equation obtained 
from (J2J is 

r + W 2 (t)r = (/ + ^)n+^xn, (6) 

P P 2 

which reduces to the planar isotropic TDHO when / = 0. Latter on, as we 

shall see, the HK method proposes a dimensional reduction from 2D to ID, 

so that we will recover the one-dimensional TDHO. 

To proceed, let us introduce the following rescaling, 

X = x/p, Y = y/p, T= f dt'/p 2 (t') . (7) 

Jo 

With the new variables, we obtain from (J5J) the rescaled stochastic equation 

^ = OxR + F(T)n, (8) 

where R = (X, Y, 0) and F(T) is defined by 

F(T) = p(t) f(t) . (9) 

Since / is a stochastic function, so is F. The statistical properties of F 
follows from (JUGjJ). The new stochastic function has zero mean, 

(F(T)) = 0, (10) 

and has correlation function 

(F(T)F(T')) = h5{T- T) , (11) 

defining a white noise. 

Equation (fTUj) is an immediate consequence of the zero mean of /. On 
the other hand, (|11J1 is demonstrated in the following way. By the definition 
of F, we have 

(F(T)F(T')) = p(t)p(t>) {fit) fit')) , (12) 



where the transformed times are 

rt dt" /■*' dt" 

/o ~pHf) ' " Jo ~pHt 



T ~ Jo 7M' ~ Jo 'oHV r )- (13) 
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From (|12j) and the correlation function in terms of the original variables, we 
obtain 

(F(T)F(T')) = hp 2 (t)5(t-t>). (14) 

This gives an expression for the correlation function of F in terms of the 
original time variable. Our objective is to show that the right-hand side of 
(fTTj) coincides with this. 
It turns out that 

using the definition of rescaled time. The right-hand side of the last equation 
can be handled with the following property of the delta function, 

^MHEfgf (16) 

for an arbitrary function (p and for (p(ti) = 0, ip'(ti) ^ 0. The sum (fTtjj) is 
over the zeros of the function ip. Applying (fTT)]) to (|To]). we get 

5{T-T) = p\t)5{t-t') ) (17) 

thus showing the equivalence between (fTTj) and (JHJ) as we desired. 

In the new variables, Equation (JHJ) is the same as that used by HK in 
the case of the time- in dependent harmonic oscillator (see Eq. (3) of 6J). 
Hence, our remaining task is to repeat the procedure by HK and map our 
conclusions to the original, non-rescaled variables. We also offer some extra 
details on the necessary calculations. 



3 Formal Solution 

We can obtain the formal solution for Equation (JHJ) as follows. Thinking the 
stochastic function F as an ordinary function, we easily obtain the linear 
invariants 

X = XcosttT + YsinttT 

1 f T 

- -^ dT' F(T') (cos fiT' + sin QT') , (18) 
V2 Jo 

Y = -XsinnT + YcosnT 

1 r T 

- -^ dT' F (T') (cos QT' -sin QT') (19) 
V2 Jo 



for the dynamical equations. In other words, both X Q and Y Q are constant 
along trajectories and can be interpreted as the initial conditions for X and 
Y respectively. Solving (J18H19|) for X and Y we obtain 



cos QT r T 

X = X cosnT-Y smQT+ ^ dT'T(T') (cos ft T' + sin ft T') 

V2 Jo 

sin ft T r T 
- — ^- / dT'FiT') (cos ft T' - sin ft T') , (20) 

V2 Jo 

cos ft T r T 
Y = y cosftT + X shift T + ^ dT'FiT') (cos ft T' - sin ft T') 

\/2 Jo 



sin ft T /" T 
+ " / dT'FiT') (cos ft T' + sin ft T') . (21) 

V2 Jo 

We can use the exact solution to expand X and Y around T = 0. Inte- 
grating (JTTJl we get the estimate F 2 (AT) ~ h/AT for small AT. Therefore, 
/ F(T)dT ~ \JhAT . Using this and expanding the exact solution, the 
result is 



X = X -QY AT + ^ dTF(T)+0((AT) 3/2 ), (22) 

V2 Jo 

1 r AT 
Y = Y + QX AT+^ dTF(T) + 0((AT) 3/2 ). (23) 

V2 Jo 

Hence we deal with an expansion in powers of s = (AT) 1 / 2 . More precisely, 

X = X + ea 1 + e 2 a 2 + 0(e 3 ), (24) 

Y = Y Q + ep 1 + e 2 [3 2 + 0(e 3 ), (25) 



where 



1 r AT 
sax = e/?i = -= / dTF{T) , (26) 

V2 Jo 

£ 2 a 2 = -ftF AT, e 2 f3 2 = nX AT. (27) 



4 Mother Field and Stochastic Quantization 

Following HK, the next step is the introduction of an auxiliary field, de- 
pending on the solution for the stochastic equation (jHJ). Let us define a 



time-independent mother field U(R) = (E/i(R), L^(R)) supposed to satisfy 
the Cauchy-Riemann conditions 

dUi dlh dUi _dU 1 

dX 3Y ' dY dX ' [ ' 

At time T, we can construct a time-dependent field Uy(R, T) which at T = 
equals U(R) and which is given in terms of the flow of the stochastic 
differential equation (|HJ), 

U r (R,T) = U(R(T)), (29) 

where R(T) evolves in time according to (JBJ). Using (J24II25J1 . we can expand 
the time-dependent field around T = to obtain 



rAT 



_ r£±± _ 

U T (R,AT) = U(R) + AT (Ox R-V)U(R)+ / dTF(T)(n- V)U(R) 

Jo 

1 rAT rAT 

+ - dT dT'F(T)F(T')(n-V) 2 U(R) 

2 Jo Jo 

+ 0((AT) 3/2 ). (30) 



where V is the gradient operator in R coordinates. 

By construction, ^^^(R, T) is a stochastic function. However, an ordinary 
function tJ(R, T) can be build after performing the last step in the HK 
method, i.e., averaging over stochastic processes, 

U(R,T) = (U T (R,T)> . (31) 

Averaging ()30|) and taking into account the statistics of F, we obtain 

§-6u, (32, 

with the time evolution operator 

H= ((OxR)-V) + |(n-V) 2 . (33) 

As shown by HK, the averaged mother field is harmonic for all times. 

Equation (}3*2j) is equivalent to a Schrodinger equation after restricting to 
the Y = line. To show this, notice that the averaging process preserves 



the harmonic property (|28|) . Using the Cauchy-Riemann properties, we can 
rewrite (J32j) as 

dUl - n(x dU2 +Y dUl ) hd2{j2 (M) 

W - ~ Q [ X -dX +Y -dx)-2VC^ (34) 

du 2 _ ( du 1 du 2 \ hd*u x 

Restricting to Y = and defining the complex field 

<KX,T)=expf-^-^J (u x {X^T)+iU 2 {X,Q,T)) , (36) 

we obtain the Schrodinger equation for a one-dimensional time-independent 
harmonic oscillator, 

with frequency Q. 

To obtain the quantization of the TDHO, consider the transformation 

^ = p- 1/2 exp(|^, (38) 

where p is a solution for the Pinney equation l{3"|). Using the Pinney equation 
and (|57)l and inverting the rescaling (J7J), we get 

. 3V> / ft 2 <9 2 w 2 (£)x 2 \ , 

'^ = (-2-^ + -y- )*■ (39) 

The spectrum of the quantum one-dimensional TDHO ()39j) can be obtained, 
for instance, solving the Pinney equation and applying the Lewis-Riesenfeld 
method (THj. 



5 Conclusion 

We have obtained the quantization of the TDHO using the HK method. Our 
approach relies on a rescaling transformation, which removes the explicit 
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time-dependence from the classical stochastic Newtonian equation. We may 
observe that, since the quantization procedure depends on the function p(t), 
which is any arbitrary solution for the Pinney equation, we have in fact an 
infinite family of quantizations. 

Since, essentially, our procedure was based on a rescaling transformation, 
it remains the question of what other classes of explicitly time-dependent 
mechanical systems are amenable to the HK method. In addition, we have 
not touched upon the measurement theory associated to HK or the physical 
origin of the noise. Finally, in the case of HK there are other possible and 
necessary extensions, like consideration of higher dimensional cases, many 
particle systems and the inclusion of general electromagnetic fields. With 
the treatment of many particle systems we can address another fundamental 
property of quantum mechanics, namely, entanglement. 
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